LBNL-61101 
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Central collisions of gold nuclei are simulated by several existing models and the central net baryon 
density p and the energy density e are extracted at successive times, for beam kinetic energies of 
5-40 GeV per nucleon. The resulting trajectories in the (p, e) phase plane are discussed from the 
perspective of experimentally exploring the expected first-order hadronization phase transition with 
the planned FAIR at GSI or in a low-energy campaign at RHIC. 



PACS numbers: 25.75.-q, 25.75.Nq. 



I. INTRODUCTION 

One of the major goals of high-energy heavy-ion re- 
search is to explore the properties of strongly interacting 
matter, particularly its phase structure [ij. The regions 
of temperature and baryon density that can be accessed 
depend on the collision energy. Thus systems with a 
very small net baryon density but rather high tempera- 
ture are formed at RHIG (satat ~ (200 GeV)^), while 
it is expected that the creation of the highest possible 
baryon densities would occur at more moderate collision 
energies (sjvtv — (6 GeV)2 0), such as those becoming 
available at the planned FAIR Q or at the low-energy 
end of RHIG. 

Our understanding of the QGD phase diagram is best 
developed at vanishing chemical potential, /is = 0, where 
lattice QGD calculations are most easily carried out. The 
most recent results indicate that the transformation from 
a low-entropy hadron resonance gas to a high-entropy 
quark-gluon plasma occurs smoothly as the temperature 
is raised, with no real phase transition being present 

On the other hand, at zero temperature most models 
predict the occurrence of a first-order phase transition 
when the density is raised , though no firm results are 
yet available for the corresponding value of the chemical 
potential, fiQ. However, if the T = transformation is 
in fact of first order, one would expect the phase bound- 
ary to extend into the region of finite temperature and 
terminate at a certain critical endpoint, {^c,Tc) @. In- 
deed, recent lattice QGD results Q suggest the presence 
of such a first-order phase transition line and an associ- 
ated critical end-point, though its precise location is not 
yet determined. 

The Facility for Antiproton and Ion Research (FAIR) 
under construction at GSI in Germany will make it pos- 
sible to create compressed baryonic matter in the labo- 
ratory, matter with a high net baryon density. The in- 
creasing interest in this physics area is underscored by 



the recent proposal for a low-energy campaign at RHIG 
aimed at the identification of the critical point Q and by 
current discussions about the feasibility of searching for 
the mixed phase at the JINR Nuclotron . 

In order to assess the prospects for using these facil- 
ities to explore the phase structure it is important to 
know what thermodynamic environments are being gen- 
erated in the bulk of the collision systems at the various 
bombarding energies available. For this purpose, we have 
employed a number of existing models to simulate cen- 
tral collisions of gold nuclei in the beam energy range 
anticipated at FAIR (5 — 40 GeV /A) and then extracted 
key information about the bulk environments generated 
in the course of a collision. 



II. THE INFORMATION EXTRACTED 

Generally, the systems involved in a high-energy nu- 
clear collision evolve rapidly in time and, furthermore, 
they are far from being uniform in space. The former 
feature prevents equilibrium from being fully established, 
while the latter feature invalidates the familiar thermo- 
dynamic relations which pertain to bulk matter. As a 
consequence, dynamical simulation models are indispens- 
able in the exploration of these processes. 

However, within a given microscopic transport model, 
it is possible to extract the characteristics of the local 
environment at any point in space and time and on this 
basis ascertain the degree of local equilibrium achieved 
and extract the corresponding local characteristics. 

Such a study was first made by Dorso et al. |lOl | who 
calculated the breakup of an initially compressed and 
heated nucleus and extracted its thermodynamic phase 
evolution. This analysis showed that the bulk of the sys- 
tem entered the spinodal region associated with the first- 
order nuclear liquid-gas phase transition. Furthermore, 
the resulting fragmentation pattern exhibited signs of fil- 
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amentation, a general characteristic feature of spinodal 
phase decomposition. 

As an instructive reference case for our present study, 
we consider central collisions of two gold nuclei and fo- 
cus on the physical conditions at the center of the system. 
Thus, in the CM frame, we consider only a small region 
around the origin, r = (0, 0, 0), and then seek to charac- 
terize the physical environment there as it evolves in the 
course of time. 

We are particularly interested in the net baryon den- 
sity p{t) = pB{t) — Psit) and the energy density e{t). 
The local stress tensor is also of interest but will not be 
examined here. Since we focus on the center of the sys- 
tem, there is no collective flow by symmetry. (While this 
is strictly true only on the average, each individual event 
might display some flow at the origin, but this possibility 
is unimportant and may be safely disregarded.) We do 
not wish to engage in a technical discussion of how these 
quantities can be extracted in the various models but 
refer the reader to the relevant literature cited for each 
particular model. We only note here that in numerical 
treatments that employ a grid in position space, such as 
fluid dynamics, these values can simply be read off at the 
appropriate lattice site (the origin, in the present case), 
while methods that represent the dynamical state of the 
system in terms of individual (test) particles must re- 
sort to an average over a suitably small test volume AV 
around the origin. [Due to the strong Lorentz contraction 
early on, the test volume must initially be sufficiently 
thin; a typical choice would be |a;|, |y|, 7cm|^| < 2fm, 
where 7cm is the Lorentz factor associated with the ini- 
tial nuclear motion in the CM frame. On the other hand, 
as the longitudinal expansion progresses and the system 
grows increasingly dilute, it may improve the sampling 
statistics to stretch the test box.] 

It is important to note that both p and e have well- 
defined values at all times. In particular, their extrac- 
tion does not require that local thermal equilibration has 
been reached. This is one advantage of considering these 
particular observables for the present study. However, 
of course, their thermodynamic relevance does depend 
on the degree of local equilibration achieved, as reflected 
principally in the isotropy of the pressure tensor. 

In each individual model, it may be possible to also 
extract local thermodynamic quantities, such as temper- 
ature T, chemical potential p, or entropy density ct, but 
although most extraction methods can be cast in suffi- 
ciently general terms to make them applicable also to 
non-equilibrium scenarios, those quantities have physical 
meaning only in equilibrium. Furthermore, importantly, 
even if equilibrium is reached, identical values of p and 
e will generally lead to different values for those thermo- 
dynamic quantities from one model to the other, due to 
their differences in mean fields and degrees of freedom. 

By contrast, the mechanical quantities p and e are in- 
herently more robust variables since they are subject to 
local conservation laws. For example, in ideal fluid dy- 
namics the conservation of four-momentum is expressed 



as d^T^^" = 0, while the conservation of baryon charge is 
expressed by the continuity equation d^j^ = 0. Since the 
various dynamical models generally abide by these basic 
conservation laws, they will have a tendency to yield sim- 
ilar results for the corresponding quantities. By asking 
about the behavior of such conserved observables we may 
therefore expect to obtain relatively robust answers. [Of 
course, for the purpose of discriminating between models 
(which is not our purpose here), it would probably be 
better to consider observables that are more sensitive to 
the specific model ingredients.] 

To further underscore the qualitative difference be- 
tween "mechanical" or "dynamical" quantities such as 
p and e and "thermodynamical" quantities such as p 
and T, we note that the above-mentioned conservation 
laws guarantee that the local energy or (baryon) charge 
density cannot change without the occurrence of a suit- 
able amount of energy or charge transport, which re- 
quires some time. By contrast, there are no such con- 
servation laws restricting the rate of change of the local 
temperature or entropy or chemical potential, which can 
change essentially instantaneously as a result of local re- 
action processes, such as ionization or chemical bonding. 
[While such phenomena may well offer useful signals of 
the hadronization phase transition (which can be viewed 
as some sort of bonding), they arc not of interest for the 
present study.] 

A further advantage of considering the variables e and 
p rather than T and p is that the equation of state, i.e. 
the pressure p{e, p) is then always a single-valued func- 
tion while this is not always the case for p(T, p) . Indeed, 
precisely when a first-order phase transition is present, 
the bulk pressure {i.e. the pressure of a spatially uniform 
system), p{T, p), is multi- valued throughout the region of 
phase coexistence. Because of this feature, if the phase 
trajectory were represented as {p{t),T{t)), it would ex- 
hibit a rather complex behavior as the expansion drives 
the system through the phase coexistence region, thus 
complicating the analysis considerably. This problem is 
not encountered in the {p{t), e{t)) representation, where 
the phase trajectory has a regular behavior throughout. 
This makes it easy, for example, to see how long time will 
be spent in the spinodal phase-coexistcncc region where 
bulk matter is mechanically unstable. 



III. PHASE DIAGRAM 

The most relevant features of the p — s phase plane 
are depicted in Fig. ^ At a given density p, the zero- 
temperature compressional energy, eT=o{p), provides a 
lower bound on the energy density e, so the accessible 
region is correspondingly limited. A useful reference is 
provided by the phase-coexistence boundaries associated 
with a recently constructed equation of state that 
has a first-order phase transition at all baryon densi- 
ties. Also shown is the corresponding boundary of a 
schematic equation of state that has a critical point at a 
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FIG. 1: [Color online] The p-e phase plane used for repre- 
senting the extracted dynamical evolution of the central en- 
vironment in Au-f Au collisions. The energetically accessible 
region is bounded from below by the zero-temperature com- 
pressional energy density er^oip)- The hadronic freezeout 
is indicated at the lower left. The phase coexistence region 
obtained in Ref. ITll| on the basis of an excluded volume is de- 
lineated by the outer contour, while the inner contour, which 
depicts a schematic boundary with a critial point, will serve 
as a reference for the phase trajectories. 



finite density, as is now generally expected. While these 
boundaries are only approximate and illustrative, they 
may serve as convenient references on the plots of the 
calculated phase trajectories. 

Fi gure Q also shows where the hadronic freeze-out oc- 
curs p^ . This representation brings out the fact that 
the freeze-out environments are quite different from those 
near the phase coexistence boundary, thus underscoring 
the importance of studying the propagation and survival 
of any proposed phase-transition signals through freeze- 
out. 



A. Isentropic expansion 

In order to establish an instructive framework for un- 
derstanding the results obtained with the various dy- 
namical models, we consider first adiabatic expansions. 
For this we use the hadronic chiral flavor-SU(3) model 
[T3| . This model is based on a chiral hadronic SU(3) La- 
grangian that incorporates and couples the complete set 
of baryons from the lowest flavor-SU(3) octet, the entire 
multiplets of scalar, pseudoscalar, vector, and axial vec- 
tor mesons, as well as baryon resonance states [TSIT3.IT5| . 
These hadrons have various types of interaction that en- 
dow them with effective masses and induce spontaneous 
chiral symmetry breaking as well as scale breaking via a 
dilaton field. The parameters of the model are fixed by 
symmetry considerations, hadronic vacuum observables, 
or nuclear matter saturation properties. The model pro- 
vides a satisfactor y d escription of both finite nuclei and 
neutron stars 0, ITgI and, furthermore, it has been 



used for fluid-dynamical studies of the space-time evolu- 
tion and HBT radii in relativistic nuclear collisions [T^ . 

With the baryon resonance couplings chosen suitably, 
the phase diagram of the SU(3) model is in qualitative 
agreement with the picture obtained from lattice results 
[l8l |. as illustrated by the corresponding phase bound- 
ary in Fig. 13 (But it is seen to differ quantitatively 
from the schematic reference boundary discussed above 
- a useful reminder of that fact that the phase boundary 
is still rather poorly understood.) Furthermore, phase 
trajectories that are consistent with the phase diagram 
of the model can be obtained by performing adiabatic 
expansions. Such expansions conserve the entropy per 
net baryon and that condition in turn yields a unique 
trajectory in the p — s phase plane, once the initial phase 
point {pi,ei) has been specified. 

A simple but rough estimate of the initial conditions 
can be obtained by assuming that the very early dy- 
namics is dominated by the interpenetration of the two 
Lorentz-contracted nuclei. Then the early baryon density 
(in the CM frame) is pi = 2^cmP{), where po w 0.15 fm~'^ 
is the normal nuclear saturation density present in the 
nuclear interior and 7cm is the Lorentz factor of the nu- 
clei in the CM frame, ^'^^^ = 1-1- Eo/2mN, where Eq 
is the beam kinetic energy per nucleon for a stationary 
target. In the CM, the energy per baryon is jcm'nT'N, 
so the energy density is = jcmruNPi = 2j^^mNPo = 
{2mN + Eo)po. The resulting phase trajectories (p, e) are 
depicted in Fig. |21 They are straight lines through the 
phase-coexistence region, while they are slightly convex 
above it and slightly concave below it. 

As we shall see, the corresponding adiabatic compres- 
sion (obtained by following these phase trajectories in the 
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FIG. 2: [Color online] The chiral hadronic SU(3) model: The 
energetically inaccessible region (grey area) and the phase- 
coexistence region (shaded area within solid boundary) are 
shown together with the corresponding quantities for the 
schematic equation of state displayed in Fig. Q Also shown 
are the phase trajectories (p, e) resulting from adiabatic ex- 
pansions starting from counterstreaming Lorentz-contracted 
nuclei at the indicated beam kinetic energies -Eo (GeV/A). 
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3-fluid: p(0,0,0,t) for Au + Au (b=0) 
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FIG. 3: [Color online] The time evolution of the net baryon 
density p{t) (top) and the total energy density e{t) (bottom) 
at the center of a head-on Au+Au collision for various bom- 
barding energies (indicated in GeV/A), in the 3-fluid model. 



opposite direction) are remarkably similar to the calcu- 
lated dynamical paths through the early non-equilibrium 
stage when counterstreaming dominates, while the subse- 
quent dynamical expansion trajectories generally exhibit 
gentler slopes. 



IV. DYNAMICAL RESULTS 

We have employed a number of different dynamical 
models in this comparitive study. Since they have been 
described in the literature already we present only brief 
characterizations here and concentrate on the resulting 
phase trajectories. 



A. Three-fluid hydrodynamics 

We first consider the three-fluid model 0] which treats 
two baryon-rich fluids originating with the incoming nu- 
clei and a baryon-free fluid created through the collisions 
among the (EoS-dependent) constituents of the first two 
fluids. The evolution of the baryon-free fluid is delayed 
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FIG. 4: [Color online] The phase trajectories {p{t),e{t)) for 
the 3-fluid collisions addressed in Fig. |3] together with the 
schematic reference phase boundary depicted in Fig. 



by a formation time r, during which it neither thermal- 
izes nor interacts with the baryon-rich fluids. After its 
formation, it starts to interact with the baryon-rich fluids 
and quickly thermalizes. With a purely hadronic equa- 
tion of state, this model was used to carry out a system- 
atic analysis of various observables at incident energies 
between few and about 160 GeV/A, and a comparison 
with results of transport models has been made as well. 

A large body of data has been well reproduced, in- 
cluding proton and pion rapidity distributions, proton 
transverse-mass spectra, A and A rapidity distributions, 
protons and pion elliptic flow (except for the proton i>2 
at 40 GcV/A), multiplicities of pions, positive kaons, 
(j) mesons, hyperons, and antihyperons, including mul- 
tistrange particles. This agreement is achieved at the 
expense of substantial enhancement of the interflow fric- 
tion, as compared to that estimated from free hadronic 
cross sections. Problems were met in reproducing the 
transverse flow ^3 1 ^-d- the directed flow requires a softer 
EoS at top AGS and SPS energies. This failure appears 
to suggest that the employed purely hadronic equation of 
state is too hard and thus leaves room for softening due 
to deconflnement. Further studies are in progress |2Cll |. 

The evolutions of the central values of p and e in head- 
on Au+Au collisions obtained with the three-fluid model 
are depicted in Fig. O for beam kinetic energies ranging 
from 5 to 40 GeV/A. At each beam energy, p{t) and 
e{t) are approximately proportional and exhibit a rapid 
growth as the two Lorentz-contracted nuclei interpene- 
trate, followed by a somewhat slower decrease reflecting 
the subsequent expansion. As the beam energy is in- 
creased, the entire history is being compressed in time. 

The separate evolutions p{t) and e{t) are then com- 
bined in Fig. 0] to yield the corresponding dynamical 
phase trajectory {p{t),e{t)). At each beam energy, the 
return path depicting the expansion lies below the early 
(outwards) path, although the two paths generally differ 
only relatively little. As was the case for the adiabatic 
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results considered above, each such "common" path is 
fairly straight and its slope increases steadily with the 
beam energy. For the lowest energy. 5 GcV/A, the tra- 
jectory just makes it to the hadronic boundary of the 
schematic phase coexistence region, while the next en- 
ergy, 10 GeV /A, already produces a turning point on the 
plasma side beyond the schematic phase coexistence re- 
gion. Thus, of the trajectories shown, this one spends the 
longest time traversing the phase coexistence region dur- 
ing the expansion phase, and the crossing time becomes 
ever shorter as the beam energy is raised. 

It should be recalled that the projection onto the (p, e) 
phase plane does not require any assumption about local 
equilibrium (which is very convenient for the purpose of 
the present study). In fact, generally, the central con- 
ditions are far from equilibrated during the early part 
of the collision. Therefore, before any thermodynami- 
cal implications could be made it would be necessary to 
carefully analyze the degree of equilibrium attained at 
any particular time of interest. 



B. Parton-hadron string dynamics 

We now consider a number of microscopic transport 
models. The first one is a recently extended version of 
the HSD model ^ called PHSD (Parton-Hadron-String- 
Dynamics) [23 |. This version includes additionally an 
early partonic phase with an equation of state from lat- 
tice QCD and quasi-particle properties for quarks, anti- 
quarks and gluons that have been obtained from fits to 
lattice results . On the hadronic side it treats explic- 
itly the familiar baryon octet and decouplet and selected 
higher resonances as well as their antiparticles. On the 
meson side it includes the pseudo-scalar and vector meson 
nonets as well as some higher meson resonances (aietc). 
Hadrons of even higher mass are treated as "strings" that 
reflect the continuum excitation spectrum of mesons or 
baryons. Since the results from the novel PHSD are very 
similar to those from HSD (without partonic phase) for 
central Au+Au collisions below about 25 GeV /A we omit 
a more detailed description of the PHSD model here. 

We recall that the HSD model has been extensively 
compared to experimental data from nucleus- nucleus re- 
actions for energies of 1 GeV /A to the top RHIC ener- 
gies {snn ~ (200 GeV)^) for hadrons made up from the 
light u, d quarks, strange hadrons as well as open 

and hidden charm |2fil |. It compares rather well with 
data (and the UrQMD model described below) for ob- 
servables such as hadron rapidity distributions |23 | but 
falls somewhat low in the A'+/7r+ ratio at top AGS or 
FAIR energies. Collective flow observables (ui (y), 1)2(1/)) 
are reproduced rather well in the SIS/ AGS energy regime 
[27II due to momentum-dependent scalar and vector self- 
energies for the baryons. However, the transverse slope 
of kaons and antikaons is underestimated for bombarding 
energies above about 5 GeV /A in central Au-|-Au colli- 
sions which has lead to the suggestion that a 'new phase 
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PHSD: p(0,0,0,t) for Au + Au (b=0) 



i i i i 1 i i i i 1 


i i i i 1 i i i i 1 i i i i 1 i i i i 
40 


/ \ 


— 30 : 




20 




10 


.... 1 .... 1 


.... 1 .... 1 .... 1 .... ■ 



4 6 8 

Elapsed time f (fm/c) 



10 



12 




Net baryon density p(t) (fm ' 



FIG. 5: [Color online] The time evolution of the central 
net baryon density p(t) (top) and the corresponding phase 
trajectory {p{t),e{t)) (bottom) at the center of a head-on 
Au-f Au coUision at various bombarding energies (indicated 
in GeV/ A), in the PHSD model, together with the schematic 
reference phase boundary depicted in Fig. Q 



of matter' or 'partonic degrees of freedom' should already 
be encountered at top AGS energies p3 |. 

The results obtained with PHSD are shown in Fig. |S1 
(Here and for the subsequent models, we do not show e(t) 
since it is approximately proportional to p{t).) The time 
evolution of the densities p{t) are remarkably similar to 
those obtained with the 3-fluid model discussed above, 
but it can be seen that PHSD leads to somewhat smaller 
compressions, particularly at higher collision energies. 

The PHSD phase trajectories are therefore also rather 
similar to those obtained with the 3-fluid model, except 
for somewhat smaller compressions and excitations at 
the highest energies. (We may also note that both mod- 
els yield a curious double-hump structure of the density 
maxima, particularly at the higher energies.) Further- 
more, at each energy, the inwards path is rather similar 
to the outwards path, though the differences are larger 
than those obtained with the three-fluid model. We also 
note that the 5 GeV /A phase trajectory turns around 
just at the hadronic border of the schematic phase coex- 
istence region (as for the three- fluid model), while the 
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10 GcN /A phase trajectory turns around at just the 
schematic plasma phase boundary of this region. 



C. Ultra-relativistic quantum molecular dynamics 

The Ultra-relativistic Quantum Molecular Dynamics 
model [2^ Isoj is a microscopic model used to simulate 
(ultra)relativistic heavy-ion collisions in the energy range 
from BEVALAC and SIS up to AGS, SPS and RHIC al- 
lowing for a consistent calculation of excitation functions. 
Its main goals arc to gain understanding of the various 
physical phenomena within a single transport model, in- 
cluding creation of dense hadronic matter at high tem- 
peratures, properties of nuclear matter, A and resonance 
matter, mesonic matter and of anti-matter, creation and 
transport of rare particles in hadronic matter, creation, 
modification and destruction of strangeness in matter, 
and emission of electromagnetic probes. 

The initial projectile and target nuclei are modeled ac- 
cording to the Fermi gas ansatz and the nucleons are 
represented by gaussian shaped density distributions. 
UrQMD includes in its collision term 55 different baryon 
species with masses up to 2.25 GeV, 32 meson species 
(including strange meson resonances) which arc sup- 
plemented by their corresponding anti-particle and all 
isospin-projected states. All these states can be produced 
in string decays, s-channel collisions or resonance decays. 
For excitations with masses higher that 2 GeV a string 
picture is used. The hadron-hadron collisions are per- 
formed stochastically like in the cascade models. The el- 
ementary cross-sections are fitted to available pp, -Kp data 
and the isospin symmetry is used whenever possible in 
order to reduce the number of individual cross-sections. 
For the interactions where no experimental data exist 
{e.g. hyperon-baryon resonance scattering), the additive 
quark model is used. 

The interactions are based on a non-relativistic 
density-dependent Skyrme-type equation of state with 
additional Yukawa and Coulomb potentials at low en- 
ergies. However, no potentials were used in the present 
calculations. For the high energy regime and baryon- 
antibaryon annihilation, UrQMD uses a string model 
similar to the Lund model [sjls^l- The strings, or the 
color tubes are first formed in the high energy hh inter- 
actions and then fragment into hadrons and new strings 
according to the Lund fragmentation procedure. For the 
ultra-high energies (top SPS energies and beyond) the 
formation of jets is also introduced into the model. 

The UrQMD results arc shown in Fig. El While the 
time evolutions of the densities p{t) and e{t) are qualita- 
tively similar to those obtained with the other models dis- 
cussed, the UrQMD compressions are somewhat higher 
than those of PHSD and more similar to the QGSM re- 
sults (see below). A likely reason for this is that neither 
UrQMD nor QGSM has any constraint on the closest 
approach between two baryons, whereas both the 3-fluid 
model and PHSD have some repulsion. 



UrQMD: p(0,0,0,t) for Au + Au (b=0) 




FIG. 6: [Color online] The time evolution of the central 
net baryon density p{t) (top) and the corresponding phase 
trajectory {p{t),e{t)) (bottom) at the center of a head-on 
Au-f Au collision at various bombarding energies (indicated in 
GeY / A), in the UrQMD model, together with the schematic 
reference phase boundary depicted in Fig. The symbols on 
the curves are separated by At = 1 fm/c. 



D. Quark-gluon string model 

The Quark-Gluon String Model (QGSM) H H HI 
Is^ incorporates partonic and hadronic degrees of free- 
dom and is based on Gribov-Reggc theory (GRT) [s^l 
accomplished by a string phenomenology of particle pro- 
duction in inelastic hadron-hadron collisions. To describe 
hadron-hadron, hadron-nucleus, and nucleus-nucleus col- 
lisions, the cascade procedure of multiple secondary in- 
teractions of hadrons was implemented. The QGSM in- 
corporates string fragmentation, formation of resonances, 
and rescattering of hadrons, but simplifies some nuclear 
effects (for example, it neglects the mean fields). 

As independent degrees of freedom the QGSM includes 
the octet and decuplet baryons, the octet and nonet 
vector and pseudoscalar mesons, and their antiparticles. 
The initial momenta and positions of nucleons inside the 
nuclei are generated in accordance with the Fermi mo- 
mentum distribution and the Woods-Saxon density dis- 
tribution, respectively. Pauli blocking of occupied final 
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states is taken into account. 

Strings in the QGSM can be produced as a result of the 
colour exchange mechanism or, like in difFractive scatter- 
ing, due to momentum transfer. The Pomeron, which is a 
pole with an intercept ap{0) > 1 in GRT, corresponds to 
the cylinder- type diagrams. The inclusive spectra in the 
QGSM have automatically the correct triple-Reggc limit 
for the Feynman variable a; — > 1. double- Rcgge limit for 
X 0, and satisfy all conservation laws. The particular 
stages of the collision model, namely (i) initialization of 
interacting projectile and target nuclei, (ii) string forma- 
tion via inelastic nucleon-nucleon (hadron-hadron) inter- 
action, (iii) string fragmentation, i.e. hadronization, and 
(iv) hadron-hadron rescattering, are solved basically by 
Monte Carlo simulation techniques. 

The results obtained with QGSM are shown in Fig. [7| 
The time evolution of the densities p{t) and e(t) are qual- 
itatively similar to those obtained with the three models 
discussed above and quantitatively closest to UrQMD, for 
the reason explained above (no repulsion). However, es- 



pecially at the lower energies, the QGSM expansion tra- 
jectories fall significantly below those of the other models. 

The model assumptions of UrQMD and QGSM are 
rather similar and the absence of hadron repulsion re- 
sults in higher maximum compressions than those pro- 
duced by the other models. However, these two models 
have significantly different expansion dynamics, probably 
due to larger number of baryonic resonances included in 
UrQMD. As a result of this, the QGSM dynamics is dom- 
inated by the propagation of nucleons and pions while 
UrQMD leads to significant resonance production. Since 
the relatively abundant QGSM pions (which are light) 
will propagate at velocities much higher than the UrQMD 
resonances (which are heavy), they will leave the inter- 
action region earlier, thus causing the energy density to 
decrease more rapidly than in UrQMD. This mechanism 
is especially important at the lower energies where heavy 
resonances play a significant role and it gradually sub- 
sides as the energy is raised, consistent with the results. 



QGSM: p(0,0,0,t) for Au + Au (b=0) 
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FIG. 7: [Color online] The QGSM evolution of the central 
net baryon density, p(t) (top), and the corresponding phase 
trajectory {p{t),e{t)) (bottom) at the center of a head-on 
Au-|-Au coUision at various bombarding energies (indicated 
in GeV/A), with the time increments At between the sym- 
bols indicated in parentheses. Also shown is the schematic 
reference phase boundary depicted in Fig. 



E. Nuclear Boltzmann equation 

The nuclear Boltzmann equation has proven to be 
quantitatively useful for the description of nuclear col- 
lisions at lower energies, up to a few GcV / A and it may 
therefore be of interest to employ it also here. For this 
purpose, we use the Boltzmann-Uhling-Uhlenbeck model 
developed by the group in GieBcn (GiBUU) [H Hi, IH 

GiBUU explicitly propagates 9 A^* and 9 A resonances 
with masses below 2 GeV as well as the S = —1 hyper- 
ons A and E and 19 hyperon resonances; the cascades 
and charmed baryons are included. The included mesons 
are: tt, 77, p, a, w, r]' , (j), rjc, J/ip, K, K, K*, K* . The 
baryon-baryon (meson-baryon) collisions below y^s = 2.6 
(2) GcV arc treated within the resonance scenario, while 
the string model is applied above. 

Thus GiBUU contains a larger set of the baryonic 
resonances than most other transport models (except- 
ing UrQMD and the Tiibingcn QMD model) and it 
consequently leads to higher pion numbers in vacuum. 
Medium corrections to the cross sections A^A^ ^ NR 
and A^A^ NNt: reduce the pion number in medium. 
The in-medium reduced cross sections are implemented 
(optionally) in GiBUU. They are computed with the 
Dirac masses from the NL2 model particular, the 

A^A^ A^A matrix element is given by the one-pion ex- 
change model, as was done in the calculations of Dmitriev 
et al. but with the vacuum A and nuclcon masses 

replaced by the Dirac values which causes a strong in- 
medium reduction of the cross section |4Cll |. 

GiBUU provides a good reproduction of nucleon col- 
lective flows ^3 as well as pion and kaon multiplicities 
mEi, at SIS energies. From AGS to the lower SPS 
energies, GiBUU overestimates pion multiplicities (with 
vacuum cross sections) but gives a reasonable description 
of the kaon multiplicities [41|, as do HSD and UrQMD. 
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GiBUU: p(0,0,0,t) for Au + Au (b=0) 
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FIG. 8: [Color online] The time evolution of the central 
net baryon density p{t) (top) and the corresponding phase 
trajectory {p{t),e{t)) (bottom) at the center of a head-on 
Au+Au collision at various bombarding energies (indicated in 
GeV/A), in the GiBUU model, together with the schematic 
phase boundary shown in Fig. Q 



The direct comparison with the HSD and UrQMD calcu- 
lations on the pion and kaon production at 2-40 GeV/A 
|4l] | has demonstrated that the model yields a some- 
what higher if "'"/tt"*" ratio due to additional meson-meson 
channels in KK production. 

GiBUU is suitable not only for nucleus-nucleus and 
hadron-nucleus collisions but also for photon-, electron-, 
and neutrino-induced reactions. This gives the possibility 
to test the same dynamical part of the model with various 
physical initial conditions. A new numerical realization 
of the model is currently being tested. 

The results presented here are based on the old ver- 
sion |3a.l39l l41|. with the calculations being done in the 
cascade mode, i.e. without a mean field and using vac- 
uum cross sections. These results are shown in Fig. |S1 
They are rather similar to those obtained with PHSD, 
the main differences being that GiBUU reaches slightly 
higher densities and the difference between the outward 
and the inward trajectory grows somewhat faster as the 
collision energy is raised. 



V. COMPARISONS 

We now compare the phase trajectories obtained 
with the different models at various collision energies 
through the anticipated FAIR range. These compar- 
isons are shown in Fig. |3 for the beam energies Eq = 
5,10,20,40 GeY/A. 

First of all, we note that by and large there is a remark- 
able degree of agreement between the results of the dif- 
ferent models. The most notable exception is the QGSM 
expansion paths which come out significantly lower than 
those of the other models, as we have already discussed 
in Sect. HVT?] 

For the subsequent discussion, in order to make a con- 
crete analysis possible, we take the adopted reference 
phase boundary at face value. But it is important to 
keep in mind that this particular boundary, though not 
inconsistent with any information we presently have, is 
likely to be quantitatively inaccurate. However, its qual- 
itative form is expected to be correct and our comments 
below are therefore expected to be robust, provided that 
appropriate adjustments are made in the specific energy 
values mentioned. 

At the lowest beam energies (under 5 GcV/A or so 
with the adopted phase boundary), the degree of com- 
pression and agitation attained does not suffice to bring 
the central part of the system into the phase coexistence 
region and such collisions are not likely to have a bearing 
on the possible existence of a phase transition. However, 
due to their relative slowness, these collisions achieve of 
high degree of local equilibrium and the data obtained 
in this range may well provide quantitative information 
on the equation of state at the corresponding moderate 
compressions. 

Above those "subthreshold" energies follows a range 
of beam energies (approximately 5—10 GcV/A for the 
adopted reference phase boundary), within which the 
highest degree of compression occurs within the region 
of phase coexistence. As the beam energy is increased 
through this range, the turning point of the phase trajec- 
tory moves across the coexistence region, starting at the 
hadronic phase coexistence boundary and ending at the 
plasma boundary. Though somewhat more violent, these 
trajectories are generally expected to still attain a high 
degree of local equilibrium. Furthermore, importantly, 
they spend the longest period of time within the phase 
coexistence region. Therefore, this energy range appears 
to be especially well suited for generating signals of the 
phase transition. 

As the collision energy is increased further, the turning 
point of the phase trajectory moves further inside the 
plasma region and, at the same time, the expansion path 
steepens. The time spent crossing the phase coexistence 
region then decreases, both in absolute terms and relative 
to the overall expansion time, so one would expect any 
phase-transition signals to gradually subside. 

Ultimately, beyond a certain critical collision energy 
(for which the expansion phase trajectory passes straight 
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FIG. 9: [Color online] The phase trajectories {p{t),e{t)) at 
the center of a head-on Au+Au collision for various bombard- 
ing energies as obtained with the indicated models, together 
with the schematic phase boundary shown in Fig. Q and the 
hadronic freeze-out line. The symbols on the UrQMD and 
QGSM curves are separated by the indicated time intervals. 



through the critical point), the phase trajectory never 
enters the coexistence region but passes entirely to the 
left of the critical point. Though interesting in its own 
right, this super-critical region of collision energy would 
not be expected to elucidate the character of the de- 
confinement phase transformation, i.e. to help determine 
whether there is in fact a first-order transition at suffi- 
ciently high baryon density. 

We also note that the adiabatic expansion results ob- 
tained with the SU(3) model correspond approximately 
to the backtracking of the very early (and mostly non- 
equilibrium) dynamics when the two Lorentz-contracted 
nuclear densities are being forced to interpenetrate. This 
feature mostly reflects the fact that the adopted initial 
values were taken to reflect such a scenario. If suitably 
modified initial conditions were chosen, for example ob- 
tained from the turning point of a dynamical trajectory, 
then the resulting expansion path would exhibit a large 
degree of resemblance with the corresponding dynamical 
trajectory. Thus one may characterize the actual dynam- 
ical expansions as being approximately adiabatic. 

Finally, we wish to emphasize that none of the dynam- 
ical models employed (except possibly PHSD) incorpo- 
rate a first-order phase transition. They would therefore 
not be suitable, in their present form, for studying actual 
dynamical consequences of a phase transition. However, 
the presence of such a phase transition is not expected 
to have an overwhelming effect on the gross dynamics, 
primarily due to the predominance of the overall expan- 
sion. [This expectation is supported by comparisons be- 
tween HSD (which does not contain a partonic phase) 
and PHSD (which does have a partonic phase) in the en- 
ergy range considered here] Therefore, it must also be 
expected that the effects of a phase transition would be 
relatively subtle and might best be studied with carefully 
designed correlation observables. 



VI. CONCLUDING REMARKS 

The present study has sought to elucidate the bulk 
conditions that may be expected to occur in nuclear col- 
lisions in the energy range where a possible first-order 
hadronization phase transition would be encountered. 
For this we have employed a number of existing dynam- 
ical models to central Au-|-Au collisions and extracted 
the time evolution of net baryon density p and the en- 
ergy density e at the center of the system where these 
are expected to achieve their largest values. The differ- 
ent models exhibit a large degree of mutual agreement 
on the behavior of the corresponding phase trajectories 
(p{t), e{t)), as was summarized in Fig. |5| even if they 
differ substantially in other regards. 

A central issue in the physics of strongly interacting 
matter is whether the hadronization phase transforma- 
tion is of first order at sufficiently high baryon density. 
The calculation of the corresponding critical point, and 
the associated phase boundary, poses a significant the- 
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oretical challenge and the question may ultimately have 
to be settled by experiment. 

The experimental investigation of this question may 
employ conceptually different strategies. One strategy 
searches for the critical point by means of special signals 
that may occur if the phase trajectory reaches its vicin- 
ity. However, the location of the critical point remains 
theoretically poorly understood and, as is well illustrated 
by Fig. 121 a small shift in its position would require a rel- 
atively large change in the collision energy of the critical 
phase trajectory that passes through it. Consequently, it 
is hard to predict what energy range would be most suit- 
able for this approach. Indeed, our present studies cannot 
rule out that the critical collision energy lies somewhere 
in the SPS range above the energies reachable by the 
planned FAIR. 

A different experimental strategy seeks direct evidence 
of the first-order transition by concentrating on signals 
that might appear as a result of the phase trajectory en- 
countering the phase-transition line. One would expect 
that such signals would best be generated if the bulk of 
the system were brought well inside the phase coexistence 
region, where a phase decomposition is favored, and kept 
there for a time sufficient to allow the development of the 
macroscopic non-uniformities associated with the phase 
decomposition. Though also associated with significant 
uncertainties, this issue can probably be assessed with 
somewhat larger confidence. Thus, considering the large 
degree of mutual agreement between the different dynam- 
ical results and taking first the adopted schematic phase 
boundary at face value, the present study suggest that 
the optimal beam energy is around 10 GeV/A, corre- 
sponding to y^J/vw ~ 2.36 GeV-|-2.36 GeV for a collider. 
We must, however, make allowance for the fact that the 
adopted schematic phase boundary is not expected to 



be quantitatively accurate. Furthermore, the dynami- 
cal models, though yielding fairly similar results, may 
all possess common inaccuracies. With a factor of two 
admitted to account for such uncertainties, the present 
study would then suggest that the optimal conditions for 
exploring the hadronization phase transition are likely 
(though not certain) to occur for beam kinetic energies 
of 5 — 20 GeV/A, corresponding to ^/snn ~ 3.6—6.4 GeV 
for a collider. 

These numbers suggest that the planned FAIR at 
GSI is well matched for such studies. Furthermore, this 
region may also be accessible at the low-energy end of 
RHIC at BNL, as well as at a possible upgraded Nu- 
clotron at JINR in Dubna. Of course, further dynamical 
studies are required before it is possible to identify the 
specific candidate signals and to assess whether they can 
indeed be expected to develop sufficiently even at the 
optimal collision energy. We hope that this study will 
provide stimulation in this regard. 
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